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ABSTRACT 

Aims. We have performed a timing and spectral analysis of the X-ray pulsar SWIFT J1626.6-5156 during a major X-ray outburst in 
order to unveil its nature and investigate its flaring activity. 

Methods. Epoch- and pulse-folding techniques were used to derive the spin period. Time-average and pulse-phase spectroscopy were 
employed to study the spectral variability in the flare and out-of-flare states and energy variations with pulse phase. Power spectra 
were obtained to investigate the periodic and aperiodic variability. 

Results. Two large flares, with a duration of ~450 seconds were observed on 24 and 25 December 2005. During the flares, the X-ray 
intensity increased by a factor of 3.5, while the peak-to-peak; pulsed amplitude increased from 45% to 70%. A third, smaller flare of 
duration ~180 s was observed on 27 December 2005. The flares seen in SWIFT 11626.6-5 156 constitute the shortest events of this kind 
ever reported in a high-mass X-ray binary. In addition to the flaring activity, strong X-ray pulsations with P spm = 15.3714 ± 0.0003 s 
characterise the X-ray emission in SWIFT 11626.6-5 156. After the major outburst, the light curve exhibits strong long-term variations 
modulated with a 45-day period. We relate this modulation to the orbital period of the system or to a harmonic. Power density spectra 
show, in addition to the harmonic components of the pulsation, strong band-limited noise with an integrated 0.01-100 Hz fractional 
rms of around 40% that increased to 64% during the flares. A weak QPO (fractional rms 4.7%) with characteristic frequency of 1 Hz 
was detected in the non-flare emission. The timing (short X-ray pulsations, long orbital period) and spectral (power-law with cut off 
energy and neutral iron line) properties of SWIFT J1626.6-5 156 are characteristic of Be/X-ray binaries. 
Conclusions. 

Key words. X-rays: binaries - stars: neutron - stars: binaries close 
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^ ■ 1 . Introduction cm 2 ). Finall y, an INTEGRAL spectrum in the range 4-60 kev 

ra/iiTrri^c <i« « *a* *a n u n ™< obtained bv Plarana et all d2006l) gave F = 0.6. £ cutoff = 6.7 keV. 

u I I ™ n T 5 ° n DecemberJ^OOS The ^ counte r p ar t is believed to be 2MASS 1 6263652- 

£ . byjhe SWIFT satellite during an X-ray outburst d Palmer et al. | 5l563Q5 fe ea et al.ll2006^ gueruela & Marco 2006), a Be star 

|2Q05|). The source showed X-ray pulsations with period ~ showi s ^ emission A chmdm observation confirmed 

15 s and strong flaring activity The spectrum was con- ^ ^ Mon rf ^ x source h consistent with the _ 

sistent with a power law with T * 3.2, giving = gested counterpart. Chandra coordinates are RA: 16h26m36.5s, 

6.1 x 10- 9 erg cm 2 s 1 in the energy range 15-100 I keV. Dec . _ 51d56m30 . 7Si with an error of 6 arcsec (Hom an, private 

Subsequent observations (see Table yj showed that SWIFT communication) 
J1626. 6-5156 is a highly va riable source with a pulse pe- 
riod of 15.37682+0.00005 s dMarkwardt & Swank! [2005b and 
pulse fraction that incr eases from ~50% between flares to 
-80% during the flares dBelloni et alj|2006l) . The X-ray spec- 
trum is well fitted by an absorbed power-law component and 

a high-energy cutoff. An iron line at 6.4 keV i s also detected In this work we present the results of a timing and spectral 
(Mar kwardt & Swank! l2505t iBelloni et all 120061) . However, the analysis of SWIFT J1626.6-5156 using RXTE data, limited to 
values of the spectral parameters differ considerably between the three observations early in the outburst where flares are ob- 
observations. [Markwardt & Swankl dlOOlh gave = 5 keV served. The details of the observation are explained in Sect. [2] 
and N H = 2 x 10 22 cm' 2 , whereas lBelloni et alj J2006I) gave and a lo g presented in Tabled In Sect. [3] the periodic and ape- 
£ cutoff = 12 keV and 7V H = 5 x 10 22 cm- 2 . ICampana et all d2006[ ) riodic time variability of SWIFT J 1626.6-5 156 is investigated, 
did not detect the iron line nor a cutoff. These authors reported Sect.0deals with the time-average and pulse-phase spectral vari- 
a harder and less absorbed spectrum (T = 0.7, 7V H = 0.95 x 10 22 abilit y- In Sect. [5]a comparison of the timing and spectral prop- 
erties of the flare and out-of-flare emission is performed and the 
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Table 1. X-ray observations of SWIFT J1626. 6-5156 in the literature and best-fit spectral parameters of an absorbed power law 
with cutoff. 



Date 


Satellite 


Flux" 


Energy range 


r 


^cut 


N b 
"n 


Fe line 


Reference 








(keV) 




(keV) 


(keV) 




18-12-2005 


SWIFT 


6.1 


15-100 


3.2 








Palmer et al. (20051 


19-12-2005 


RXTE 


4.1 


2-40 


0.9 


5 


2 


6.4 


Markwardt & Swank (2005) 


19-12-2005 


RXTE 


3.6-2.4 


3-25 


0.9-1.3 


12 


4-5.5 


6.4 


Belloni et al. (2006) 


12-01-2006 


SWIFT 


1.2 


0.5-10 


0.67 




0.94 




Campana et al. (2006) 


21-01-2006 


INTEGRAL 


1.6 


4-40 


0.6 


6.7 






Tarana et al. (2006) 



a: xlO 9 erg cm 2 s 1 
b: xlO 22 cm" 2 
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Fig. 1. Light curve of Obs. 1 showing the X-ray pulsations and 
the increase of the pulse fraction during the flare event. Bin size 
is 1 s. The dashed lines mark the flare interval used in the analy- 
sis (see text). 



2. Observations and data analysis 

SWIFT J1626.6-5156 was observed with the Proportional 
Counter Array (PCA) onboard the Rossi X-ray Timing Explorer 
(RXTE) starting from December 19, 2005, one day after the 
discovery by Swift. Here we concentrate on the three observa- 
tions made between December 24-27, 2005 (MJD 53728.385- 
53731.295), in which the source displayed flaring activity. The 
data correspond to RXTE proposal P91094. Four PCUs were 
on during the entire duration of the three observation intervals. 
PCU1 was off all the time. 

Light curves and spectra were extracted using version v.5.3.1 
of the RXTE FTOOLS package. The following screening crite- 
ria were applied to the data prior extraction: good time inter- 
vals were defined when the pointing of the satellite was stable 
(< 0.02"), the elevation above 10° and far away from the South 
Atlantic Anomaly. 

The log of the observations is given in Table|2] Each observa- 
tion contains one flare event, although the event that took place 
in Obs. 3 was significantly smaller and shorter. The duration of 
the flares varied from -180 to 480 s. The flare in Obs. 2 was de- 
tected at the end of the observation interval. Thus, the duration 
of post-flare interval of Obs. 2 is short (~ 100 s). The pre-flare 
interval of Obs. 3 contains a data gap of about 2260 s due to the 
passage through the SAA anomaly. 




PCU2 2-15keVc/S 



Fig. 2. Light curves and hardness-intensity diagrams of the three 
observations. The Y-axis was left the same in the three panels to 
allow easy comparison. The entire time interval is used with a 
bin size of 64 s. The dashed lines indicate the limits of the pre-, 
post- and flare intervals selected for the analysis. 

3. Timing analysis 

3.1. Light curves and hardness ratios 

Strong pulsations are detected at the 15.37s period from all ob- 
servations. Figure Q] shows the light curve of Obs. 1 at high res- 
olution (At = 1 s). The top panel shows the first 1.2 ks of data, 
where the pulsation is evident. The bottom panel shows the sec- 
ond part of the observation, which includes a strong flare. Its 
shape is very similar to that of the other strong flare (from Obs. 
2): a major increase in the pulsed fraction, clearly accompanied 
by a more moderate increase of the DC level, followed by a de- 
pression of DC level and pulsed fraction and a slow recovery. 
Fig.|2]shows the three flare events with corresponding hardness- 
intensity diagrams using a coarser time bin (At = 64 s). The 
hardness ratio was defined as the ratio between the count rate in 
the 9-15 keV band to the count rate in the 2-9 keV band. The flare 
in Obs. 2 is very similar to that in Obs. 1 , although it appears later 
in the observation and the post-flare region is not observed. The 
flare in Obs. 3 is much weaker. As the count rate increases, the 
hardness ratio increases until the count rate reaches about 350 c 
s at which point it flattens. The dashed lines in Fig. [T]and|2] 
indicate our definition of pre-flare, flare and post-flare intervals 
for each observation. 

As mentioned, during the flares, the pulsed fraction increases 
considerably; the minima of the pulsation cycle increase by a 
factor of 2, while the maxima increase by a factor of 5. After 



Table 2. Summary of the observations. 
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Obs. 


RXTE ID 


MJD 


MJD 




Pre-flare 






Flare 






Post-flare 




interval 




start 


end 


'x 


Duration 


PF 


/ ° 

i peak 


Duration 


PF 


'x 


Duration 


PF 


1 


91094-02-02-01 


53728.385 


53728.425 


995 


1800 


44% 


3555 


450 


72% 


750 


1142 


31% 


2 


91094-02-02-02 


53729.168 


53729.211 


988 


2900 


46% 


3880 


480 


70% 


600 


108 


36% 


3 


91094-02-02-04 


53731.228 


53731.295 


680 


2057 


38% 


1040 


180 


51% 


645 


1236 


37% 



a: Background-subtracted count rate for 4 PCU in the energy range 2-15 keV for a time bin of 16 s. 
Durations are in seconds. 



the flares, both the average rate and the pulsed fraction go below 
their pre-flare values, and start a slow recovery. In Obs. 1 the flux 
had not completely recovered the pre-flare level 1200 seconds 
later, when the observation ended. Given the fact that the pre- 
flare level of Obs. 2 is similar to that of Obs. 1, it is likely that 
the source intensity fully recovered to the pre-flare level. 

In order to investigate differences in the emission properties 
between the flare and quiescent states, we isolated the time inter- 
vals corresponding to the flare events. For each observation, we 
defined three intervals: pre-flare, flare and post-flare. The exact 
definition of these intervals can be seen in Fig. [2] 

The average count rate (peak intensity in the case of the 
flare interval), duration and pulse fraction of the available in- 
tervals corresponding to the three "states" (pre-flare, flare and 
post-flare) for each observation are given in Table [2] 

3.2. Spin period determination 

In order to determine the spin period of SWIFT J1626. 6-5156 
we obtained a 0. 125-s binned barycentric corrected light curve of 
the non-flaring emission for the three observation intervals in the 
energy range 2- 1 5 ke V (from PC A BJ2ms.8BDJ5_Q mode) . An 
FFT applied to the light curve revealed a coherent modulation at 
0.065 Hz, which correspond to a pulse period of 15.4 s. Several 
of its harmonics are clearly seen. In order to accurately measure 
the pulse period, a pulse-folding analysis was performed. The 
light curve was divided into blocks of continuous stretches of 
data points with a duration of 230 s. Each block was folded mod- 
ulo the initial period found from the FFT. The resulting pulse 
profile obtained for the first block was used as a template. The 
difference between the actual pulse period at the epoch of obser- 
vation and the period used to fold the data was determined by 
cross-correlating the pulse profiles corresponding to each block 
and the template. The phase shifts were fitted to a linear func- 
tion, whose slope represents the correction in frequency to the 
actual period. A new template was derived for this corrected pe- 
riod, and the process repeated. The derived pulse periods were 
15.3718+0.0003 s, 15.371 1+0.0002 s and 15.3724±0.0007 s for 
the three observations, respectively. 

3.3. Pulse profile and pulse fraction 

The average pulse profiles at three different energy ranges, ob- 
tained by folding the corresponding light curves onto the best-fit 
spin period, are shown in Fig. [3] The pulse profile of the flare 
is somehow narrower than that of the out-of-flare but for both 
there is no strong dependence on energy. The pulse fraction was 
calculated as PF =(/ mflX -/ m ,„)/(/ mflA - +/„„•„). 

In order to perform a more in-depth analysis of the individ- 
ual pulse profiles, we extracted a light curve from the full PCA 
energy range with a bin size of P sp i n /20, divided it in intervals 




Pulse phase Pulse phase 



Fig. 3. Comparison of the pulse profile as a function of energy 
during the flare (right) and non-flaring (left) part of the light 
curve. Data correspond to Obs 1 . 

20-points long (centered approximately on the pulse maximum) 
and with these constructed a time-phase image, with time on the 
X-axis, pulse phase on the Y-axis and count rate on the Z axis. 
The image on a gray scale color map is shown in Fig.|H together 
with the corresponding light curve (top panel). From this fig- 
ure, one can see that the out-of-flare pulse peaks at around phase 
0.5, while the flare itself is characterized by a maximum around 
phase 0.6. Moreover, at T~500 one can see a small rate increase 
corresponding to a slightly higher phase, indicating that a very 
small flare is probably present there. 

In order to study in more detail the pulse changes during the 
flare, we accumulated an out-of-flare template by simply adding 
the first 100 cycles. Comparing this template with the single 
pulses during the flare, it appears as if the flare could be charac- 
terized by the appearance of one or two additional components 
at fixed phases. 

In order to assess this numerically, we ran the following 
procedure. Each cycle was fitted with the out-of-flare template 
with a variable overall multiplicative factor. Not surprisingly, 
outside the flare and post-flare region (and with the exception 
of three cycles around cycle 35), this model gives a decent fit 
with normalization equal to 1 . By decent fit we mean a fit with a 
X 2 r = ^ 2 /20 <15: this is not formally acceptable, but no system- 
atic residuals are seen and given the high statistics it can not be 
expected that a single simple model could fit all single cycles. 
Overall, out of flare the^j; values are around 8. 

For the cycles with;^ > 15, we identify the largest positive 
deviation from the template fitting (limited to the phase range 
0.2-0.8) and repeat the fit with the addition of a Gaussian com- 
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Fig. 4. Top panel: light curve of Obs. 1. Bottom panel: corre- 
sponding time-phase image, obtained aligning all pulses as a 
function of time. Darker areas indicate higher count rate. 
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Fig. 5. Values of xj — X 2 /20 obtained with the fitting procedure 
described in the text. The different curves correspond to simple 
template (dotted line), template plus one Gaussian (dashed line) 
and template plus two Gaussians (full line). The horizontal line 
marks the level x 2 r = 15. 



Fig. 6. Best fit parameters as a function of time. On all three pan- 
els, the gray curve shows the PCA count rate evolution for ref- 
erence. Top panel: multiplicative factor for the out-of-flare tem- 
plate. The dashed line shows the level of 1. Middle panel: phase 
of the centroid of the Gaussians (first Gaussian: filled circles; 
second Gaussian: empty circles). Bottom panel: normalization 
of the Gaussians (first Gaussian: filled circles; second Gaussian: 
empty circles). The dashed line corresponds to zero normaliza- 
tion. 

ues for five of the seven parameters are shown in Fig. [6] The 
remaining two parameters, the width of the two Gaussians, are 
clustered around 0.1-0.2 and 0.01-0.05 in phase for the first and 
second Gaussian respectively. 

Based on these results, we can identify three intervals: 

- Before the flare: the pulse shape is consistent with a constant 
shape, corresponding to what we used as template. 

- During the flare: the template shows a higher normalization 
factor (up to 3), while an additional broad Gaussian cen- 
tered at phase 0.6 appears. Occasionally, a second narrower 
Gaussian is necessary for a fit. 

- During the post-flare depression: the template normalization 
is lower than unity and a negative Gaussian with a very con- 
stant normalization is apparent. 

- After the post-flare depression: the template normalization 
slowly recovers to unity, no Gaussians are needed. 



ponent centered at the position of that maximum deviation, with 
normalization equal to the maximum deviation and with start- 
ing width 0. 1 . With this addition, only a few cycles around the 
maximum of the flare retain a^ > 15. 

In order to try to find a better fit for these cycles, we added 
an additional Gaussian component. The procedure is the same as 
before. The new Gaussian is added only if the fit with a single 
Gaussian results in a^ > 15. The fits are considerably better 
during the flare, although there are still a few points with a xl > 
15. The^ values obtained with this procedure are shown in Fig. 
® 

The final fits have a number of free parameters between one 
(the template normalization factor for the template-only fits) and 
seven (template normalization plus three parameters for each of 
the two Gaussians) for the most complex fits. The best fit val- 



3.4. Power spectrum and noise components 

In order to obtain the power spectra, the Obs. 2 light curve was 
binned with bin size 2~ 9 s, which is the maximum resolution pro- 
vided by the BJ2msSBJ)35jQ data mode. Then the light curve 
was divided into segments of 512 s for the non-flare part and 
128 s for the flare. An FFT was calculated for each segment (6 
for the non-flare and 4 for the flare). The final power spectra re- 
sulted after averaging all the individual power spectra and rebin- 
ning in frequency. Figure [7] shows the power spectra and noise 
components of the pre flare and flare light curves. In addition to 
the harmonic components of the pulsations, strong band-limited 
noise is detected both in the quiescence (i.e. non-flaring emis- 
sion) and during the flare, with an integrated 0.01-100 Hz frac- 
tional rms of around 40% and 64%, respectively. The change in 
the pulse profile between pre-flare and flare (as seen in Fig. [3]) 
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Fig. 7. Power spectrum of SWIFT J1626.6-5156. Coherent pul- 
sations are clearly seen at ~0.065 Hz together with some higher 
harmonics. The quiescence emission shows a weak QPO com- 
ponents at 1 Hz. The lines indicate the Lorentzian components 
used to fit the power spectra (see text) 
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Fig. 8. 5-d binned ASM light curve of SWIFT J1626.6-5156 for 
the period MJD 50135-54340. The start time of the PC A obser- 
vations is indicated by an arrow. 



can also be expected from the change in the relative strength of 
the harmonics of the pulsation. 

The power spectra were fitted with Lorentzians. Ignoring 
the X-ray pulsation and its harmonics, a total of 5 Lorentzians 
were required to fit the non-flare power spectrum, while the flare 
power spectrum needed 3 Lorentzians. 

Three Lorentzians are common to the two data sets. These 
are zero-centred Lorentziants that fit the continuum. The power 
spectra of the non-flare part require two more Lorentzians. One 
is centred at the same frequency as the pulsation and can be 
interpreted as the coupling betwe en the periodic and aperiodic 
variability dLazzati & Stellal[l997l) . This coupling can be recog- 
nised by the broadening of the wings of the narrow peaks of 
the periodic modulation. The lack of this component in the flare 
power spectrum maybe due to the fact that the aperiodic noise is 
stronger (see the higher continuum short-ward of the pulsation) 
during the flare. The non-flare power spectrum presents a QPO 
at around v = 1 Hz with Q (v/FWHM) of 4 and fractional rms 
amplitude of 4.7%. 



Fig. 9. After the outburst that led to its discovery as an X-ray 
source, SWIFT J1626. 6-5156 exhibited a periodic modulation 
with period ~45 days. On the right the hardness ratio-intensity 
diagram is shown. 



3.5. Longer-term variability 

Figure [8] shows the long-term variability of SWIFT J1626.6- 
5156 and the transient nature of its X-ray emission. The source 
spends most of the time in a quiescent state (the average count 
rate is *0.06 ASM cs 4 ) and only occasionally goes into out- 
burst. Note that significances above 3 sigma are only seen in the 
5-day average ASM bins that cover the 2005 December outburst. 

We also collected the average PCU2 count rates for each of 
the 295 RXTE observations, thus resulting in an unevenly spaced 
time series (see Fig.[9j. The onset of the PCA observations coin- 
cided with the outburst peak (Fig. [8). A first visual inspection of 
the resulting lightcurve shows two main properties: i) an initial 
fast flux decay lasting for about 160 days since the outburst on- 
set, followed by ii) a series of apparently periodic "bumps". In 
order to carry out a first estimate of the two components we fitted 
the first 160 day points with an exponential component and the 
second part of the lightcurve with a sinusoid. For the first com- 
ponent we found that r of the exponential is ~33 days. For the 
second component we obtained a period of 45+1 days, though 
the reduced^ 2 is only about 2. 

We note that bumps are rather asymmetric, suggesting either 
that the sinusoid model is a rough estimate of the raw data and/or 
the 45-d period is a harmonic of the true period. This interpre- 
tation is also corroborated by the very high asymmetric profile 
of the latest four peaks, with the odd and even bumps differing 
each other. To test this hypothesis we add an additional sinu- 
soidal component to the fit with a period forced to be twice to 
that of the first component. The inclusion of this component has 
a probability of ~99.Q% to be significant. It is evident that only 
future optical spectroscopic studies of the companion star and/or 
pulse arrival time analysis of the X-ray photons can provide an 
unambiguous detrmination of the orbital period. The first bump 
after the exponential flux decay phase is a reminiscence of what 
is routinely observed in man y Be/X-ray binaries as they enter 
in an outburst phase (see e.g. Campan a et alJll999h : a rapid rise 
in flux, a flat topped maximum and a relatively slow decay. At 
variance with the above interpretation is the second bump which 
falls at the minimum of the expected modulation. 

Figure [9] also shows the variation of the hardness of the 
source, defined as the ratio between the count rates in the en- 
ergy bands 6-15 keV to 2-6 keV, as a function of intensity. As 
the outburst decayed the source became softer . 

It is worth mentioning that the same flare mechanism that 
occurs during high-intensity states, i.e., near the outburst peak, 
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Table 3. Results of the fits to the power spectra. Only the non- 
coherent components are shown. 



Lorentzian 
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Fig. 10. Energy spectra of the flare and out-of-the-flare parts for 
Obs 1 . The lower panel shows the ratio between the two spectra. 



still works at lower mass accretion rates. The variability seen 
at the end of the outburst decay, around MJD 53800-53900 (see 
Fig-UK is caused by the same type of flaring activity that we have 
reported in previous sections. The flares show the same factor 3- 
4 increase in flux and also an increase in the pulse fraction. 

4. Spectral analysis 

4.1. Time-averaged spectroscopy 

In order to investigate potential spectral variability between the 
non-flare emission and the emission during the flare, we ex- 
tracted energy spectra corresponding to the flare and a stretch 
of non-flare emission. Fig. [TOl shows the spectra for Obs.l and 
Table 14. ll summarises the results of the spectral fits for the three 
observations. A systematic error of 1% was added in quadrature 
to the statistical one. An absorbed power law with a high-energy 
cut off and an iron line at ~ 6.4 keV provided a good fit to the 
spectral continua. No spectral variability is seen within the un- 
certainties, although the flare spectrum in Obs. 2 is somehow 
harder than the preflare one. The pre-flare and flare spectra ba- 
sically differ in that the pre-flare spectra require the inclusion of 
an edge at ~9 keV to obtain an acceptable fit. The probability 
that the improvement of the fit occurs by chance (F-test) by the 
addition of this component is <10~ 5 . 




Pulse phase 



Fig. 11. Spectral parameters as a function of pulse phase for the 
preflare state of Obs. 2. 

4.2. Pulse-phase spectroscopy 

As mentioned above, the shape of the pulse profile of SWIFT 
J1626. 6-5156 hardly changes with energy. A more detailed way 
to investigate how the X-ray spectrum changes as a function of 
the spin of the magnetized neutron star is by performing pulse- 
phase spectroscopy. In order to investigate the pulse phase de- 
pendence, high time and energy resolution are needed simulta- 
neously. Unfortunately, the only data mode available in our ob- 
servations that provides detailed time and energy information is 
BJ2ms_8BJ)_35_Q, which provide high time resolution but mod- 
est energy resolution, namely, eight energy channels in the en- 
ergy range 2.6-15 keV. 

The pulse profile of Obs. 2 was divided into eight equally 
spaced bins to provide eight 2.6-15 keV background-subtracted 
spectra. Each spectrum was fitted with a power law plus iron 
line and edge. Due to the modest energy resolution, the iron 
line parameters were not well constrained. The iron line and 
edge energy were fixed to 6.4 keV and 8.8 keV, respectively, 
the value obtained from the time-average spectrum. As can be 
seen in Fig. QT| despite that the flux changed by a factor 2, no 
significant changes were seen in any of the spectral parameters, 
neither in the continuum (power-law photon index) nor in the 
line and edge (width, T e d g e)- This result confirms the absence of 
variability of the pulse profiles with energy. 

5. Discussion 

We have performed a timing and spectral analysis of the newly 
discovered X-ray pulsar SWIFT J1626. 6-5156. We have con- 
centrated on the observations made at the early phase of the re- 
ported outburst (see Fig. [8] and |9]l in which the source displays 
flaring activity. Each X-ray light curve was divided into subin- 
tervals corresponding to the out-of-flare and flare parts of the 
observations. A separate analysis was carried out in each subin- 
terval. 

5. 1 . Out-of-flare analysis 

The out-of-flare parts of the light curves were used to derive 
the spin period because these parts do not present high ampli- 
tude flux changes and because of their longer duration, allowing 
enough number of cycles to be included. No significant varia- 
tions in the spin period were seen during the observations. Thus 
we derived a weighted mean spin period of 15.37 14±0. 0003 s. 

The pulse profile of the out-of-flare part of the X-ray light 
curve, obtained by considering 100 cycles folded on the spin 
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Table 4. Spectral fits results for the three observations. Uncertainties are 90% confidence. Flux corresponds to the energy range 
3-30 keV. 



Obs. 1 Obs. 2 Obs. 3 



p ar ameter s 


pre -flare 


flare 


pre- flare 


flare 


pre- flare 


flare 


Flux (erg cm" 2 s" 1 ) 


3.5 x 10" 9 


8.4 x 10~ 9 


3.5 x 10~ 9 


6.3 x 10~ 9 


3.1 x 10~ 9 


3.8 x 10~ tJ 


N H (xlCr 1 cm l ) 


9.4* 


9.4* 


9.4* 


9.4* 


9.4* 


9.4* 


r 


„ _ _ _i_n HQ 
"•'-'-0.05 


0.6±0.1 


0.8+0.1 


0.5+0.1 


0.82+0.06 


0.8±0.1 


norm. 


11+001 


22+0 03 


12+0 01 


1 5+0 02 


12+0 01 


1 3+0 02 


E cul (keV) 


4.7±0.3 


4.8±0.4 


5 +01 


4.5±0.4 


5.0±0.3 


5.0±0.4 


Efcu (keV) 


Q 7+0.8 
^•'-0.4 


8.6±0.6 


10.3+0.5 


o n+0.6 
S -2-0.3 


10.0 + °* 


10 6 +L4 


E Fe (keV) 


6.2±0.1 


6.5±0.2 


6.3±0.1 


6.4+0.1 


6.3+0.1 


6.4+0.1 


o-Fe (keV) 


0.4* 


0.4* 


0.4* 


0.4* 


0.4* 


0.4* 


norm. (xl0~ 3 ) 


6±1 


8+3 


5 2 +1 


9±2 


5±1 


7±1 


EW(Fe) (eV) 


195+30 


105±30 


180+30 


160±20 


195±35 


235±80 


£ edge (keV) 


S./_ .4 




8.9+0.5 




8.9+0.5 




^"edge 


0.07±0.03 




0.07±0.02 




0.06±0.02 






0.98(46) 


0.85(48) 


0.99(46) 


0.91(48) 


0.82(48) 


1.15(48) 



*: fixed 



period, and a variable overall multiplicative factor fits well the 
individual pulses with normalization equal to unity. In the case 
of the post-flare interval, where both the overall flux level and 
the pulse amplitude are depressed, the best fits are obtained with 
the addition of a negative Gaussian component. 

The average non-flare X-ray flux in the 3-30 keV band is 
3.5 x 10~ 9 erg s _1 cirT 2 for Obs. 1 and Obs. 2 and slightly lower 
in Obs. 3, 2.9 x 10~ 9 erg s cm 2 . There is no clear evidence 
for variability of the spectral parameters with time nor between 
the flare and out-of-flare spectra in any of the three observations. 
The only difference is the presence of an absorption edge com- 
ponent in the out-of-flare spectrum that it is not required in the 
flare spectrum. 

5.2. Flare analysis 

The flare event consists of a sudden increase of the X-ray flux 
by a factor of four and is accompanied by an increase of the 
pulse fraction and a hardening of the spectrum. The average 3- 
30 keV flux, taking into account the entire duration of each flare, 
decreased from 8.5 x 10~ 9 erg s _1 cirT 2 in Obs. 1 to 6.3 x 10~ 9 
erg s _1 cirT 2 in Obs. 2 and further to 3.8 x 10~ 9 erg s _1 cm -2 for 
Obs. 3. 

The template fitting showed that the flare corresponds to the 
appearance of an additional component at a different phase than 
the "normal" pulse shape outside flares, in addition to a bright- 
ening of the out-of-flare template. The average pulse profile ob- 
tained before the flare does not generally fits well the individual 
pulses during the flare. The addition of a Gaussian component 
centered around phase 0.6 (shifted from the peak at phase 0.5 
of the template shape) reduces largely (see Fig. [5]) the amount 
of bad fits. Still, there are some pulses that do not fit well to 
an average profile plus a Gaussian. The addition of a second 
Gaussian improved further the fits during the flare. In the post- 
flare phase, when the template normalization goes below unity, 
this additional component disappears and is replaced by a nega- 
tive Gaussian component centered at phase 0.5. As the template 
slowly recovers to pre-flare levels, the negative component also 
disappears. 

This behavior is difficult to understand. The flare appears 
to be made of a general increase in the pulse fraction plus the 
appearance of a peaked component at phase 0.6. The relatively 
small difference in average energy spectrum between non-flare 



and flare indicates that the spectrum of the additional compo- 
nent is not very different from that of the template emission. The 
post-flare negative contribution shows that after the flare the ad- 
ditional Gaussian disappears and is replaced by a depression of 
the template maximum, together with a decrease of the pulse 
fraction. 

5.3. The nature of the source 

The X-ray spectral components of SWIFT J1626. 6-5156, 
namely, absorbed power-law component plus high-energy cutoff 
and iron line at 6.4 keV are characteristics of accreting X-ray 
pulsa rs dWhite et alj 119831; iReig & Rochel 1 19991 ICoburn et all 
|2002|) . Although a few low-mass X-ray pulsars are known (e.g. 
Her X-l) the vast majority of pulsating neutron star are found 
in high-mass X-ray binaries. Massive X-ray binaries are classi- 
fied according to the luminosity class of the optical component 
into supergiant X-ray binaries and Be/X-ray binaries. The former 
tend to be persistent sources, have short (P 01 -b < 10 d) orbital 
periods and long spin periods (P sp i n > 100 s), while Be/X-ray 
binaries are usually transients, have wider orbits (P or b > 20 d) 
and contain fast rotating neutron stars (P sp i n < 100 s). On long 
timescales Be/X-ray binaries display two types of X-ray activ- 
ity: regular, orbitally modulated outbursts normally peaking at 
or close to periastron with X-ray flux increases of about one or- 
der of magnitude with respect to the pre-outburst state, reaching 
peak luminosities L x < 10 37 erg s _1 (type I) and bright (X-ray 
flux 10 3 - 10 4 times that at quiescence) uncorrected with orbital 
phase and that last for several orbital periods (type II). 

The long-term variability of SWIFT J1626. 6-5 156 shown in 
Fig. [8] and [9] is reminiscent of Be/X-ray binaries (Reig 2007b . 
especially, the detection of smaller periodic outbursts (type I) 
following a major one (type II). This behaviour has been see n 
in the Be/X-ra y binaries 4U 01 15+6 3 dNegueruela et al.lll998l) 
KS 19 47+300 ilGallowav et al]l2004T) and EXO 2030+375 (iReigl 
2008). The periodicity is associated with the orbital period of 
the system. The X-ray emission increases as a result of accre- 
tion from the circumstellar disk of the Be star during periastron 
passage. Therefore, it seems natural to assume that the small out- 
bursts in SWIFT J1626. 6-5156 are also separated by the orbital 
period. Nevertheless, given the asymmetric profiles of the out- 
bursts (note the last four shown in Fig. [9j, we cannot exclude 
that the true orbital period is twice the modulation. In fact, the 
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optical co unterpart to SWIFT J 1626.6- 5 156 exhibits strong Ha 
emission. Negueruela& Marco (2006) reported an Ha equiva- 
lent width of ~ -40 A. Such e quivalent w idth implies an orbital 
period 100 < P olb < 200 days (Reig 2007). Although the Ha line 
may also appear in emission in supergiant binaries, the equiva- 
lent width is normally below ~1 A in these systems. 

We can then conclude that all the observational data available 
of SWIFT J1626. 6-5156 both from the vicinity of the compact 
object (relatively short spin period, X-ray spectral parameters, 
transient X-ray emission) and the optical companion (strong Ha 
emission, relatively long orbital period) suggest the presence of 
a neutron star orbiting a Be star companion. 

If the nature of SWIFT J1626.6-5156 as a Be/X-ray binary 
is finally confirmed, then this is the first time that such super- 
fast flaring activity is reported for this type of systems. Up to 
now the fastest X-ray variability seen in Be/X-ray binaries was 
that associated with the rotation of the neutron star (X-ray pul- 
sations). Flaring activity is not uncommon in massive X-ray 
binaries habouring supergiant companions (SGXR) but absent, 
with only one exception (EXO 2030+375), in Be/X-ray bina- 
ries. During an X-ray outburst of the Be/X-ray binary EXO 
2030+375 detected by EXOSAT in 1985 October, a series of 
six flares that recurred quasi-periodically every 3.96 hr were ob- 
served. The duration of these flares varied between 1.3-2.3 hr 
and showed significant variability during the decay, consisting 
o f intensity quasi-pri odic oscillations with periods of 900-1220 
s (Parmaretal 1ll989h . 

The flares of EXO 2030+375 share similar properties with 
those seen in LMC X-4 flares, including the ~3-4 hr flare recur- 
rence time, the -1-2 hr flare duration, the high fluences of the 
flares (> 10 40 ergs), the rapid rise with gra dual decay, and th e 
absence of any significant spectral change dMoon et al.l l2003). 
These similarities suggest a common origin. One possible mech- 
anism could be Rayleigh-Taylor insta bilities near th e magneto- 
spheric boundary of the neutron star ( Appar aoll991l) . 

Flaring activity has also been reported in two wind-fed su- 
pergiant X -ray binary. 4U 19 07+09 exhibits flares on time scales 
of 1-2 hr dFritz et al.ll2006b . These flares are locked to the or- 
bital motion as they o ccur twice per orbit (lin't Zand et al.l fl998; 
Muk eriee et al.l200lb . Vela X-l has been ob served to flare above 
its persistent level o n several occasions (lLaurent etal.lll995t 
iRrivonos et al.l l2003). The duration of the flares was 5-6 hr. 

However, the duration of the flares in the systems mentioned 
above is longer, typically by one order of magnitude, than the 
flares observed in SWIFT J1626. 6-5156. It is their short du- 
ration and the fact that they occur in a system with a main- 
sequence companion what makes the event in SWIFT J1626.6- 
5156 unique. 

Fast and short-lived flares have been also reported for the 
newly suggested class of supergiant X -ray binaries, known as 
supergiant fast X- ray transients, SFXT (Neguer uela et al.l i2006; 
ISmith et al.ll2006b . As in the classical supergiant binaries these 
spora dic outbursts typically last for few hours (ISguera et al.l 
2006). They differ in that SFXTs are transient sources and in 
the complex structure of the outbursts, consisting of one or a 
few short flares. However, at least two SFXTs have shown flares 
with a duration compara ble to that of SW IFT J1626. 6-5156, 
namel y XTE J 1739-302 dSakano et al.l l2002b and IGR J 17544- 



namel y Alfc J 17 M-Ml ( Sakano et 
2619 (Gonzalez-Rie stra et al.ll2004l) . 



The physical explanation of the fast outbursts is not clear. 
Accretion from a smooth homogeneous stellar wind is ruled 
out. Instead, porous (clumpy) winds are invoked. Fast outbursts 
would be due to the accretion of large, dense clumps on the neu- 



tron star (IWalter & Zurita-Herasll2007T: lNegueruelaetail l2008). 
Recently, Sido li et ail (120071) has proposed the presence of a sec- 
ond wind component in the form of an equatorial disk around the 
supergiant donor to explain the outbursts seen in the SFXT IGR 
Jl 1215-5952, the only SFXT that exhibits periodic outbursts. 
These outbursts last for about 10 days and exhibit flaring activ- 
ity near the peak. They are reminiscent of the type I outbursts so 
common among Be/X-ray binaries and most likely have a com- 
plete different orig in to that of the fast, short (~ 1 hr) flares. Even 
ISidoli et ail (120071) suggest that the ~l-hr duration flares seen at 
the peak of the longer outburst in IGR J 1 1215-5952 can be as- 
cribed to the clumpy nature of the wind in the equatorial disk. 

Th e asymmetry of ma ss outflows in Be stars is well estab- 
lished dWaters et al.lll988l) . The structure of the stellar winds 
in Be stars contains two components: at higher latitudes, mass 
is lost through the high-velocity, low-density wind (typical of 
early-type stars); in the equatorial regions a slower and denser 
wind operates and ultimately forms the disk. The interaction be- 
tween the neutron star and the equatorial disk gives rise to en- 
hanced accretion and increase of the X-ray flux. 

The short-lived flares in SWIFT J1626.6-5156 can then be 
explained by the interaction between the accreting compact ob- 
ject and a clumpy stellar wind. Whether this stellar wind is the 
polar or equatorial component is not clear and will probably de- 
pend on the geometry of the orbit. If the equatorial plane and 
orbital plane are not coplanar and the orbit is relatively wide (so 
that the compact object does not directly crosses the disk) then 
the neutron star will be exposed to the polar wind. The determi- 
nation of the orbital parameters of SWIFT J1626. 6-5156 may 
provide the key to understand the physical mechanism behind 
the flares. 

Given the short duration of the flares in SWIFT J1626.6- 
5156 it is tempting to compare the SWIFT J1626. 6-5156 prop- 
erties with those of neutron star low-mass binaries, in particu- 
lar those of the "bursting pulsar" and the "Rapid Burster" since 
these systems also show fast and short episodes of increase X- 
ray emission. Instead of the term "flares" the word "burst" is 
used in this context. Both systems are transient, and the accretor 
in both is a neutron star. Both show postburst "dips" followed by 
a relatively long recovery interval preceding a further burst. This 
behavior can be understood in terms of depletion of the reser- 
voir followed by a fill-in time. Though the burst duration is very 
long in the case of SWIFT J1626. 6-5156 it is also interesting 
to note that burst durations up to 700 s have been detected from 
the Rapid Burster from type-II bursts. The longer bursts of Rapid 
Burster (comparable to those observed in SWIFT J 1626.6-5 156) 
have a recovery time of ~lhr, clearly not covered by the rela- 
tively short (lhr) RXTE pointings. 

Type II bursts likely originate by spasmodic release of grav- 
itational potential energy, which is almost certainly due to some 
not yet understood accretion disk instability. In support of this 
similarity of bursts is the fact that the spectra of the burst and 
the persistent emission are nearly identical, suggesting a similar 
(resca led) mechanism for both dBriggs et al.lll996t iLewin et all 
fl996T) . 

6. Conclusions 

The available observational data indicate that SWIFT J1626.6- 
5156 is a member of the class of massive X-ray binaries known 
as Be/X. The relatively short spin period, long orbital period, 
spectral components, transient emission and strong Ha emis- 
sion are typical of this type of systems. However, the short-lived 
flares make SWIFT J1626. 6-5156 unique in its class. 
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The time-average energy spectrum can be represented by an 
absorbed power law (T sa 0.7), modified at high energy by a cut- 
off (Eta ~ 4.8 keV, E foli * 9 keV. An Fe line at -6.4 keV as also 
present. The out-of-flare power spectra show a QPO at ~ 1 Hz, 
which is not present in the flare power spectra. 

The out-of-flare pulses are consistent with an average pro- 
file. During the flare (and three early cycles around cycle 35), 
the cycle can be decomposed into the enhanced template of the 
out-of-flare pulse, plus a main "Gaussian" component starting 
at phase 0.5 and drifting monotonically to phase 0.6. Its width 
is about 0.15-0.3 in phase. During the post flare, the curve is 
consistent with a de-enhanced template and a negative Gaussian 
component of about the same width as the previous one and a 
phase around 0.5. 

Although the actual triggering mechanism of the flares is 
not known, we speculate that it can be associated to a non- 
homogeneous, i.e. clumpy stellar wind. 
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